From the Cooper problem to canted supersolids in Bose-Fermi mixtures 
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We calculate the phase diagram of the Bose-Fermi Hubbard model on the 3d cubic lattice at 
fermionic half filling and bosonic unit filling by means of single-site dynamical mean-field theory 
(DMFT). For fast bosons, this is equivalent to the Cooper problem in which the bosons can induce 
s-wave pairing between the fermions. We also find miscible superfluid and canted supersolid phases 
depending on the interspecies coupling strength. In contrast, slow bosons favor fermionic charge 
density wave structures for attractive fermionic interactions. These competing instabilities lead to 
a rich phase diagram within reach of cold gas experiments. 
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Interactions between bosons and fermions play a cru- 
cial role in various physics contexts. Examples in- 
clude the atomic nucleus, quarks exchanging gluons via 
the strong force, electrons dressed by lattice vibrations 
forming polarons, conventional superconductors where 
phonons induce an attraction between the electrons at 
the Fermi energy, and the phase separation between "^He 
and "^He mixtures. Beyond mean-field, these systems 
are notoriosly difficult to describe. Cold atom experi- 
ments can be used to simulate this physics, thanks to 
the experimental control over the coupling strength be- 
tween fermions and bosons, effectively performing quan- 
tum simulation of superconductors. 

The first experiments investigated the influence of 
fermions on the bosonic Mott insulator, and found 
that the bosonic visibility always decreases when adding 
fermions attractively interacting with the bosons [l|, 12] ■ 
This has been explained by self-trapping [2|-|5|, cor- 
rections to higher bands |j, |6|, or by adiabatic heat- 
ing [lll7|,|8|. At weaker inter-species interactions, symme- 
try between repulsion and attraction was found [5| . In a 
dynamics experiment the strength of the potential terms 
has been measured with astonishing precision [9|. How- 
ever, many more exo tic p hases such as supersolids llOll 
and pair supcrfluids 11 1 have been predicted |12l4l9| . 
though not yet realized in experiment. Such may be- 
come possible though thanks to the recent discovery of 
multiple Feshbach resonances between ^^Na and "'"K at 
MIT 0. 

In this Letter, we revisit the Cooper problem of con- 
ventional superconductors in a cold atom setup, that is 
we study the conditions under which bosons induce s- 
wave pairing between spin-1/2 fermions [21ri23| . We will 
see that a bosonic condensate leads to a strong static en- 
hancement of s-wave pairing. Our formalism also allows 
us to explore physics in the strong Bose-Fermi coupling 
regime as well as bosons that are slow compared to the 
Fermi velocity. In such cases, instabilities favoring den- 



sity waves compete against pairing, leading to a rich and 
unexpected phase diagram. 

Our model consists of spinless bosons and spin-1/2 
fermions on a cubic lattice with Hamiltonian 

('U> {i'j) i<^ « 

i i icr 

where h\ and hi (cj^ and da) are the bosonic (fermionic) 
creation and annihilation operators at site i with spin 
CT and ri° {n\) denote the corresponding number oper- 
ator. Particles can hop between neighboring sites via 
the hopping amplitude ib(f) ^nd the particle number is 
adjusted through the chemical potential /ib(f)- The parti- 
cles can interact via an onsite interaction, where C/bb- % 
and Uhi denotes the boson-boson, fermion-fermion and 
boson-fermion interaction, respectively. We will work at 
unit filling for the bosons and half filling for the fermions, 
in which case the sign of t/bf is irrelevant. This model is 
a direct extension of the previous cold atom experiments 
with spin-polarized fermions. We restrict the discussion 
to the case where the spin-up and spin-down fermions 
interact equally strongly with the bosons. 

To numerically study the above model we use DMFT 
where the full many body problem is mapped onto a self- 
consistent determination of an impurity model. In the 
Nambu notation the kinetic impurity action for sublattice 
s is given by 

ST = -\ f^dTdr'{h\{T) - *t)A,^,(r - r')(b,(r') - *,) 

-Z<*1, / dThs{T)~ / dTdT'ct(T)Af,,(T-T')c,(T') 
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where Ab(f) is the matrix hybridization function of the 
bosons (fermions) and the corresponding creation and de- 
struction operators are given by b|(r) = (&|(r), ^^(t)) 



and c|(r) = (cl (r), C4^.s(r)), and s ~ A,B denote 

the two sublattices. $1^ = (b''').^ = {(f>^s,4>-s) is 
the time independent condensate order parameter of the 
bosons determined selfconsistently on the other sublat- 
tice as denoted by the subscript — s. For a cubic lat- 
tice, the coordination number is z = 2d = 6. The 
fermionic hybridization function is determined by the fol- 
lowing form of the inverse lattice Green function [c^ (r) = 



G{. \fc,iw„) 



C - Sa -Sa -Efc 

. Efe -Ss -c + sj 



(with C = iu!n + fi,€k = 2tf ^ ■ -^ cos(A:j), and stan- 
dard notation for the normal and anomalous selfenergies 
on the respective sublattices) such that (charge) density 
wave ordering and s-wave pairing are allowed, and can 
occur independently of each other. The nature of the 
density-density coupling between bosons and fermions 
implies that a density wave ordering for fermions im- 
mediately creates density wave ordering for the bosons, 
and vice versa. The possible symmetry breakings in the 
spin sector are expected to be the same as in the pure 
fermionic model. The (local) potential energy terms are 
absorbed in the potential part of the impurity action 
'S'pot = JodTHiocir). 

As impurity solver we use a continuous-time Monte 
Carlo method based on an expansion of the partition 
function in powers of the impurity-bath hybridization 



^b(f) 



and the condensate order parameter $. 
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impurity solvers. This method allows for 



method is a direct extension of the fermionic 



the first time to study Bose-Fermi mixtures within the 
fuU DMFT formalism { ^, see however Refs.fii, [26 1 
regarding the broken symmetry in the action). An illus- 
tration of a possible Monte Carlo configuration is shown 
in Fig. 2.^ Details of the algorithm will be presented else- 
where 29|. 

At half filling the pure fermionic system (C/bf = 0) 
exhibits particle-hole symmetry: The superfluid phase 
transition on the attractive side Ug < is mirror reflected 
around C/ff — into a anti-ferromagnetic transition on the 
repulsive side. Us > 0, as is shown in Fig. [2] (although 
both have SU(2) character, we already use the terminol- 
ogy appropriate for C/bf 7^ 0) . The DMFT results interpo- 
late between the Weiss mean- field result Tmf = 6i^/|C/ff| 
valid at strong coupling and the T-matrix/BCS result at 
weak coupling [28|. We first study how the superfluid 
and anti-ferromagnetic phase transition are affected by 
the presence of strongly condensed bosons with a speed 
of sound exceeding the Fermi velocity (referred to as fast 
bosons), and focus on the s-wave pairing transition. The 
bosons can then be treated in the Bogoliubov approx- 
imation I23II and the effective interaction between the 
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FIG. 1. (Color online) Illustration of a typical Monte Carlo 
configuration. The full (empty) circles denote creation (an- 
nihilation) operators in the imaginary time interval [0, /3) 
for bosons (black), and spin-up (blue) and spin-down (red) 
fermions. a) The local contribution to the weight of the op- 
erator sequence is determined by the length of the segments 
and the overlap between segments of different particles (seg- 
ments mark time intervals in which a particle resides on the 
impurity), b) A possible configuration of bosonic hybridiza- 
tion functions and source fields with density rib = at imagi- 
nary time T — 0. c) All possible combinations to connect the 
fermionic creation and annihilation operators. 



fermions is given by 



Us 
Us 



C^bfXo(k,w) 



(1) 



C/^f2nb(ztb 



^l) 



w2 - (ztb + eb)((zfb + Ek) + 2nbC/bb) 



with xo(k, tj) the density-density response function. 
With a strong condensate, the zero temperature expres- 
sion can be used since the Bose condensation tempera- 
ture is much higher than the BCS temperature. When 
the bosonic sound velocity Sh = (2nb^bbib)^'^ is much 
higher than the Fermi velocity, retardation effects are 
negligible [23[ and the induced interaction is always at- 
tractive on the Fermi sphere. The induced interaction is 

rind/'i^^ _ _Ubt_„_n^\ — t/fat 1 

C/bb 



then U'^^^ik) 



^ci(k) 



(7bb i+e{^~i:Ui 



s(fcja)) 



(with ^ = yJt}j/2nhUhh the healing length), and an on- 
site effective interaction U^ = Us — jj^ X^k "^i C^) ^^ 
found. The effective hopping follows from a mean- field 
decoupling of the nearest neighbor interaction and is 
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FIG. 2. (Color online) S-wave superfluid (left) and antifer- 
romagnetic (f/bf = and |f/bf|/tf = 8 on the right) phase 
transition of the Bose-Fermi Hubbard model on the 3d cu- 
bic lattice at filling Tib = 1 and n^ = 714, = 1/2 and with 
Uhh/tf = 20 and ib/if ~ 1 for different boson-fermion in- 
teractions Ubf- The DMFT results interpolate between the 
Weiss mean-field result and the T-matrix/BCS result ('mod- 
ified HF') (see text). Inset: Critical temperature for pairing 
for non- interacting fermions {Uff =0) as a function of the 
boson-fermion interaction Uht- The transition temperature is 
exponentially suppressed at low Ua for all Ubt (not shown). 
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This leads to the phase diagram shown in Fig. [5] where 
for small \Us\ s-wave pairing is enhanced by stronger 
boson-fermion interactions and anti-ferromagnetism is 
suppressed. Pairing can hence occur for non-interacting 
and repulsive pure fermions. The inset of Fig. [2] shows 
that the transition temperature in the purely induced 
case {Us = 0) can be of the same order as for an attrac- 



tive fermionic system without bosons [23|. This holds 
even for values of rihU^h far outside the Bogoliubov 
regime. For stronger interspecies interactions than the 
ones shown, phase separation occurs |10l . ISCl |31 | which 
prevents a further increase of Tc. 

The shape of the phase boundary for large C/bf in Fig. [2] 
looks surprisingly similar to the phase diagram of the 
purely fermionic system. On the basis of the perturba- 
tive arguments given above, we look for effective inter- 
actions t/g^ and effective hoppings tf^ of the respective 



Uff - c[U^i/Ui,b and tf - if - c'^U^i/Uy 



t 

off' 



bb 



forms Ug^ 

with c'l and C2 fitting constants. In Fig. |3] we see that 
all transition lines can be collapsed onto each other; i.e., 
that in the presence of a fast bosons self-consistent first 
order contributions suffice to explain the physics, even 
far outside the perturbative regime. An analysis of the 
quasi- particle weight [28|, Zqp = (1 — ImE(iwo)/wo)^^ 
(with Wo = ttT the lowest Matsubara frequency measured 
from the Fermi level) in systems where symmetry break- 
ing was disabled confirmed this picture (not shown) : Zero 
quasi-particle weight corresponds to the Mott insulator 
on the repulsive side {Uff > for C/bf = 0) and the molcc- 



FIG. 3. (Color online) The phase diagrams of Fig. [2] can be 
collapsed onto the phase diagram of a pure fermionic model 
with renormalized hoppings tf — tf — OiUbt/Uhb and on- 
site repulsions Ug = Uff — c'lUbt/Uhb, with ci and C2 fitting 
constants. Error bars are of the order of the symbol size and 
omitted for clarity. 



ular density wave on the attractive side. Collapse of the 
curves with different Uht is observed provided the on-site 
repulsions, hoppings and Zqp factors are rescaled. 

The density-density correlation function in Eq. ([T]) 
changes dramatically in the absence of a condensate. It 
may change sign when uj cannot be set to zero thereby 
suppressing pairing. This motivates us to numerically in- 
vestigate the dependence of the phase transition on the 
bosonic hopping ty,, shown in Fig.|3}3 for strong interac- 
tions Uff/tf = —10. We see that the system undergoes a 
sharp first order transition around th/tf w 0.75 between a 
fermionic superfluid (corresponding to spin singlets in the 
fermionic spin sector) and a (molecular) charge density 
wave (corresponding to Neel ordering in the fermionic 
spin sector). The bosons remain strongly condensed at 
this point (no ~ 0.6), but pick up charge density wave 
order. The transition temperature varies remarkably lit- 
tle over the different phases, reflecting the underlying 
SU{2)xSU{2) symmetry of the pure fermionic model. At 
very low hoppings (ib/if < 0.2), the bosons become insu- 
lating and are very ineffective in influencing the fermions. 
The fermions can undergo a simultaneous a pairing and 
molecular charge order transition, which couples back to 
the bosons and generates bosonic charge order. We also 
observed that, except in the close vicinity of a bosonic 
supcrfluid-insulator phase transition, bosonic static mean 
field approximation provides quantitatively correct re- 
sults in our DMFT scheme, which may be useful for fu- 
ture cluster extensions of this work. 

We repeated this calculation for different values of U^f 
for a temperature T/tf = 0.2 close to the ground state re- 
sulting in the phase diagram in the (f/bf, tb) plane, shown 
in Fig.|3]4_. For large values of t/bf we find the same phases 
as in Fig.|3J\: a double superfluid, a CDW with a bosonic 



superfluid, and a CDW with a fermionic superfluid. 

However, for rather low values of U^ and sufficiently 
large bosonic hoppings wc find a supersolid phase, in 
which bosons and fermions have both types of orderings. 
In this supersolid, the gaps for pairing and charge order 
are not equal; this supersolid is a realization of the canted 
supersolids put forward in Refs [32|, |33| . The RG study 
of Ref. [3J| finds that a d-wave superfluid develops for 
certain parameters in this regime, which may compete 
with the supersolid. However, seeing such a phase is not 
possible with single site DMFT. We expect d-wave only 
to be feasible for low values of Ue and C/bf while for large 
values of Ug and U^f the supersolid is most likely sta- 
ble. The transition temperature of the supersolid phase 
for \Uh{\ = 2 and ib = t[ is Tc ~ 0.48tf, rendering an 
experimental observation with cold gases realistic. This 
is the same transition temperature as for a supersolid in 
a bosonic model on a triangular lattice [35|, and 50% 
higher than the one of an anti-ferromagnct in the 3d 
Hubbard model [36[- One example of a mixture with 
promising scattering properties for the supersolid phase 
is ^Li-^Li 37|. For low values of Us the structure of the 
phase diagram is identical to the one shown in Fig. |4]A., 
from which we conclude that the BCS-BEC crossover is 
not a driving force for the Bose-Fermi Hubbard model at 
half filling. 

In conclusion, we developed a single-site DMFT for- 
malism for the Bose-Fermi-Hubbard model allowing for 
s-wave pairing and charge density wave ordering. We 
computed changes to the pure fermionic phase diagram 
at fermionic half filling induced by the commensurate 
bosons, focusing on attractive Ug- While fast bosons fa- 
vor s-wave pairing, slow bosons favor charge density or- 
der. These different type of instabilities compete, leading 
to some unexpected phases such as a canted supersolid 
and the CDW-|-SFb phase shown in the phase diagram 
of Fig. H 
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FIG. 4. (Color online) Panel A: Low temperature phase di- 
agram of the Bose-Fermi Hubbard model on the 3d cubic 
lattice in the {th,Uht) plane. The parameters are f/bb/if = 
20, f/ff/if = -10,T/ff = 0.2, n^ = n; = 1/2, nb = 1. The 
labels are: charge density wave (CDW), superfluid fermions 
(SFf), superfluid bosons (SFb) and canted supersolid (SS = 
CDW+SFf+SFb). The CDW-hSFt to SFb+f and CDW+SFb 
to CDW-f SFf are flrst order, the other ones are second or- 
der. Panel B: Critical temperature for s-wave pairing as a 
function of the bosonic hopping amplitude ib for |[/bf |/ff ~ 8 
(indicated by the dashed line in Fig. [31A_). The bosonic con- 
densate no is shown for a temperature corresponding to the 
maximum of the pairing and charge density wave ordering 
temperatures. 
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